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SYSTEM FOR METHOD FOR FINE ACQUISITION OF 
A SPREAD SPECTRUM SIGNAL 

FIELD OF THE INVENTION 

This method relates, in general, to CDMA (Code Division 
Multiple Access) spread spectrum receivers, and more 
specifically, to fast acquisition of a spread spectrum signal, 
such as is performed by a fast acquisition GPS (Global 
Positioning System) receiver. 

BACKGROUND OF THE INVENTION 

Spread spectrum communication in its basic form is a method 
of taking a data signal that is used to modulate a sinusoidal 
carrier and then spreading its bandwidth to a much larger value, 
e.g. in a global positioning system (GPS) application by 
multiplying a single -frequency carrier by a high-rate binary (- 
1,1) pseudo-random noise (PRN) code sequence that is known to 
GPS users. Thus, the signal that is transmitted includes a data 
component, a PRN component, and a (sinusoidal) carrier 
component . 

At the receiver, a synchronized replica of the transmitted 
PRN code is required to de-spread the data sequence. Initial 
synchronization, called acquisition, is followed by fine 
synchronization, which is called tracking. 

The present invention relates to acquisition. Acquisition 
is the process by which the replica PRN code is synchronized (to 
within a small timing offset) with the code conveyed by the 
received signal either for the first time or after losing a 
previously acquired signal, and also by which the carrier 
frequency of the received signal is determined. Thus, to 
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acquire a signal, an acquisition system must accurately 
determine any frequency- shifting of the received signal from the 
transmitted frequency in order to accurately wipeoff (remove) 
the carrier signal. Frequency-shifting can be caused by 
relative motion of the transmitter and receiver (Doppler- 
shifting) as well as by clock inaccuracies (so that a 
transmitter and receiver sometimes do not agree on what is in 
fact the same frequency) . The carrier frequency-shifting 
results in a modulation of a code component after carrier wipe- 
off in the receiver. Thus, in acquiring a signal, it is also 
necessary that the replica code sequence be not only time- 
aligned with the received code sequence, but also modulated to 
compensate for the frequency- shifting so as to fully eliminate 
the PRN sequence and leave behind only the data conveyed by the 
received signal. The acquisition process is therefore a two- 
dimensional search, a search both in code phase and in 
frequency. 

For low level signals additional frequency refining is 
sometimes necessary for closing tracking loops. Because of the 
two dimensional search needed for initial acquisition, a useful 
strategy is to first search with as coarse a frequency grid as 
possible in order to accelerate the search. As long as the 
replica code is aligned with the incoming signal in the code 
dimension one can refine the carrier frequency using. only a one- 
dimensional but more fine resolution search in frequency. 

A bottleneck here is the data modulation, i.e. the 
modulation of the signal by the data. In the initial 
acquisition stage, one strategy that is used is to take the 
frequency information from a coherent stage. The coherence 
length increase is limited due to 50Hz navigation message and it 
reduces the gain of the system beyond 2 0msec integration. Thus, 



-2- 



for maximum gain, the resolution of coherent processing based 
methods is limited to 50Hz. To go beyond this range one must 
either compensate for data modulation or sacrifice gain. 

In weak signal conditions, a fine acquisition stage is 
required for tracking initialization. The frequency resolution 
of initial acquisition schemes is limited, due to the complexity 
of the two-dimensional search (i.e. for code phase and for any 
shift in the carrier frequency) and due to the data rate of 50Hz 
(i.e. the navigation data rate for a GPS receiver). To obtain 
fine resolution, the data modulation must be wiped off, which 
according to the prior art has been accomplished either by bit 
detection, which is not always possible in weak signal 
conditions, or by a search over all possible data modulations of 
the signal as sampled in a set of samples, which imposes a 
significant computational burden. Another technique has been to 
use so-called zero-padding of a short segment of a signal so as 
to end up with a segment of a desired length, suitable for 
spectral analysis using a discrete Fourier transform (DFT) with 
a frequency spacing determined by the length of the zero-padded 
signal segment. Zero-padding with fine frequency spacing, 
without more however, does not yield fine frequency resolution 
in a noisy environment. 

What is therefore needed, especially in weak signal 
conditions where sacrificing gain is sometimes not practical, is 
a computationally efficient way to more finely resolve the 
carrier frequency search than is possible without taking into 
account the data modulation. 

SUMMARY OF THE INVENTION 

Accordingly, the present invention provides a method, 
corresponding apparatus, and corresponding system for acquiring 
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a received spread spectrum signal, the received signal having a 
carrier component at a carrier frequency, a code component 
having a code period, and a data component, the acquiring 
including matching the phase of a replica of the code component 
to the phase of the received code component and also determining 
any shift in the carrier frequency away from a transmitted 
carrier frequency, the method including the steps of: performing 
a first acquisition of the received signal so as to provide an 
approximately estimated carrier frequency and a phase of the 
replica and also so as to provide a code-wiped and an 
approximately carrier-wiped signal; and performing a second 
acquisition of the approximately carrier-wiped signal, the 
second acquisition including a substep of data wipe-off 
involving a squaring or similar operation on a signal derived 
from the approximately carrier-wiped signal; thereby providing a 
correction to the approximately estimated carrier frequency, a 
correction that accounts for the carrier frequency shift 
remaining after the first acquisition. 

In a further aspect of the invention, the step of 
performing a second acquisition of the approximately carrier- 
wiped signal includes the substeps of: mixing the approximately 
carrier-wiped signal with a plurality of sinusoids in turn, each 
sinusoid at a different frequency in a range of frequencies 
serving as trial corrections to the approximately estimated 
carrier frequency determined in the step of performing a first 
acquisition, so as to provide a mixed signal; performing a first 
coherent processing of the mixed signal so as to provide a 
carrier-amplified signal; performing a data wipeoff using the 
carrier-amplified signal so as to provide a data-wiped signal, 
the data wipeoff including a squaring or similar operation on 
the carrier-amplified signal, thus retaining in the data-wiped 
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signal, information about the carrier frequency shift remaining 
after the first acquisition; performing a second coherent 
processing of the data-wiped signal, the second coherent 
processing for providing a further correlated and filtered 
signal; and detecting the best value to use for the correction 
to the approximately estimated carrier frequency based on the 
further correlated and filtered signal for each different trial 
frequency correction. ^ 

In another, further aspect of the invention, the step of 
performing a second acquisition of the approximately carrier- 
wiped signal includes the substeps of: performing a first 
coherent processing of the approximately carrier-wiped signal so 
as to provide a carrier-amplified signal; performing a data 
wipeoff using the carrier-amplified signal so as to provide a 
data-wiped signal, the data wipeoff including a squaring or 
similar operation on the carrier-amplified signal, thus 
retaining in the data-wiped signal, information about the 
carrier frequency shift remaining after the first acquisition; 
and performing a mixing and second coherent processing, the 
mixing being performed on the data-wiped signal using a 
plurality of sinusoids in turn, each sinusoid at a different 
frequency in a range of frequencies serving as trial corrections 
to the approximately estimated carrier frequency determined in 
the step of performing a first acquisition, to provide a mixed 
and data-wiped signal, and the second coherent processing being 
performed on the mixed and data-wiped signal, the second 
coherent processing for providing a further correlated and 
filtered signal. 

In yet another, further aspect of the invention, the step 
of performing a second acquisition of the approximately carrier- 
wiped signal includes the substeps of: performing a first 



coherent processing of the approximately carrier-wiped signal so 
as to provide a carrier-amplified signal; performing a data 
wipeoff using the carrier-amplified signal so as to provide a 
data-wiped signal, the data wipeoff including a squaring or 
similar operation on the carrier-amplified signal, thus 
retaining in the data-wiped signal, information about the 
carrier frequency shift remaining after the first acquisition; 
and performing a cross correlation of the data-wiped signal, the 
cross correlation including estimating the phase and then the 
correction to the approximately estimated carrier frequency. 

In yet still even another, further aspect of the invention, 
the step of performing a second acquisition of the approximately 
carrier-wiped signal includes the substeps of: performing a 
common coherent processing of the approximately carrier-wiped 
signal so as to provide a first-processed; mixing the first- 
processed approximately carrier-wiped signal with a plurality of 
sinusoids in turn, each sinusoid at a different frequency in a 
range of frequencies serving as trial corrections to the 
approximately estimated carrier frequency determined in the step 
of performing a first acquisition, so as to provide a mixed 
signal; performing a first coherent channel processing of the 
mixed signal so as to provide a carrier-amplified signal; 
performing a data wipeoff using the carrier-amplified signal so 
as to provide a data-wiped signal, the data wipeoff including a 
squaring or similar operation on the carrier-amplified signal, 
thus retaining in the data-wiped signal, information about the; 
performing a second coherent channel processing on the data- 
wiped signal, the second coherent processing for providing a 
further correlated and filtered signal; and detecting the best 
value to use for the correction to the approximately estimated 
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carrier frequency based on the further correlated and filtered 
signal for each different trial frequency correction. 

The invention also comprehends a method including only a 
single acquisition, an acquisition including a data wipe-off 
involving a squaring or similar operation on a signal having a 
data component and a frequency component, i.e. the invention is 
of use in acquiring any kind of signal having a carrier 
component and a data component, regardless of whether a first, 
typically relatively coarse acquisition is performed. 

Thus, according to the present invention, to achieve better 
than 50Hz resolution, the data modulation of a CDMA signal is 
compensated for by complex squaring of the signal (or a similar 
operation) . In the case of a binary phase shift keying (BPSK) 
data signal d(t)~±\ with frequency shift modulation oxp(jwt) (i.e. 
w is the shift in the frequency of the carrier) , complex 
squaring doubles the frequency and removes the data modulation, 
essentially according to, 

s\t) = d(t)eJ Wt d(t)eJ' Wt =e 2jWt . 

Therefore, the frequency of the resulting signal is twice the 
angular frequency w sought to be determined. So by taking the 
square of the data modulated signal (not the square of the 
magnitude of the data modulated signal, but rather the ordinary 
algebraic square), the data modulation can be wiped off. 
However, such an operation is sometimes not, by itself, enough 
for robust acquisition to better than 20dB-Hz (signal-to-noise 
(S/N) ratio multiplied by the frequency band at the point where 
the S/N is taken) , although it usually is enough in case of a 
strong signal. In general, though, the squaring technique must 
be combined with properly incorporated coherent accumulation 
(i.e. accumulation of terms taking into account both their 
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magnitude and sign or phase) and filtering stages to achieve 
fine resolution better than 1Hz with only a moderate 
computational burden . 

BRIEF DESCRIPTION OF THE DRAWINGS 

The above and other objects, features and advantages of the 
invention will become apparent from a consideration of the 
subsequent detailed description presented in connection with 
accompanying drawings, in which: 

Fig. 1A is a schematic/ block diagram of a portion of a 
receiver in which, according to the present invention, the an 
initial acquisition module is followed by a fine acquisition 
module for providing a fine frequency shift correction to 
refine the carrier frequency determined by the initial 
acquisition module; 

Fig. IB is a schematic/ block diagram of a the fine 
acquisition module of Fig. 1A, showing that it includes a 
correlator followed by a fine acquisition kernel; 

Fig. 2A is a schematic/ block diagram of a first embodiment 
of a fine acquisition kernel according to the invention, an 
embodiment in which the output of the correlator is 
immediately demodulated based on a guess of the fine 
frequency shift correction, followed by a first coherent 
processing module; 

Fig. 2B is a schematic/ block diagram of the same first 
embodiment of a fine acquisition kernel as illustrated in 
Fig. 2A, but specifically for an implementation in which 
processing is performed in parallel for each channel, i.e. 
for each trial (guess) fine frequency shift correction, and 



-8- 



including a low pass filter and downsampler as an 
implementation of the first coherent processing module; 

Fig. 2C is a schematic/ block diagram of the same first 
embodiment of a fine acquisition kernel as illustrated in 
Fig. 2A, and specifically for an implementation in which 
processing is performed in parallel for each channel, as in 
Fig. 2B, but instead of using a low pass filter and 
downsampler as an implementation of the first coherent 
processing module, an integrate and dump module is used, 
followed by a filter; 

Fig. 2D is another implementation of the first embodiment 
of a fine acquisition kernel, the embodiment indicated 
generally in Fig. 2A, an implementation in which zero padding 
is used; 

Fig. 3A is a schematic/ block diagram of a second 
embodiment of a fine acquisition kernel according to the 
present invention, one in which the output of the correlator 
is not demodulated (based on a guess of the fine frequency 
shift correction) until after the first coherent processing; 

Fig. 3B is a schematic/ block diagram of an implementation 
of the second embodiment, an implementation using a discrete 
Fourier transform in which the output of the correlator is 
not demodulated (based on a guess of the fine frequency shift 
correction) until after the first coherent processing; 

Fig. 3C is a schematic/ block diagram of another 
implementation of the second embodiment, an implementation 
using a cross correlation, instead of a mixing and second 
coherent integration; and 

Fig. 4 is a 3rd embodiment of the present invention, one 
similar to the second embodiment as indicated in Figs. 3A-C, 
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with a common first coherent processing module, but one in 
which the bandwidth of the signal after correlation is 
reduced by filtering the signal. 

BEST MODE FOR CARRYING OUT THE INVENTION 

The present invention will now be described in connection 
with a receiver for receiving a direct sequence (DS) spread 
spectrum signal transmitted by a satellite vehicle. A global 
positioning system (GPS) signal transmitted by a satellite 
vehicle is an example of a DS spread spectrum signal. It should 
be understood however, that the present invention, a system for 
acquiring a spread spectrum signal, is not limited to acquiring 
a spread spectrum signal in any particular context, such as GPS, 
nor is it limited to acquiring a particular kind of spread 
spectrum signal. As will be clear from what follows, a system 
according to the present invention for acquiring a spread 
spectrum signal can be of use in any application in which a 
spread spectrum signal is to be acquired, as long as the spread 
spectrum signal includes a code component and a carrier 
component . 

Notation and context 

In what follows, it is assumed that the carrier is wiped- 
off the received signal, and we are processing the digitized 
signal after analog to digital conversion. The initially 
periodic signal is corrupted by noise and shifted in frequency 
due to the Doppler effect and other causes (including 
differences in the rates of the transmitter and receiver 
clocks) , and is modeled as 

x(n) = s(n) + x noise (n) 

where the received signal at the receiver is 
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,j(W D n+<p) 



s(n) = s 0 (n-n c )e 



= s 0 (n-n c )e 



in which 



s 0 (n) = A-d{floor(n/N dp )}r(nmodN cp ) . 



Here, s(n) is the received signal at the receiver, x noise (n) is the 
noise component of the received signal, A is a scaling factor, 
s 0 (n) is the signal with zero code phase and Doppler frequency or 
other frequency shift offsets, n c is the unknown code phase of 
the received signal (with respect to the replica code) measured 
in for example samples, N dp is the duration of a data bit (in 

samples) , N cp is the code period (in samples) , f D is the 

unknown frequency shift appearing because of relative motion of 
the receiver and the transmitter and because of other causes, 
W D is the corresponding angular frequency normalized to the 
sampling frequency, e {-l,l} is the sample value obtained in 
sampling a data bit, and r(w)e{-l,l} is a sample of the replica 
code. The variable w is used to represent the ordinary angular 
frequency, related to the ordinary frequency / by the equation 
w-2jrf . The disclosure also uses the variable W t as above, to 
represent the digital signal processing analog of the ordinary 
angular frequency. The variable w' is a normalized angular 
frequency and is dimensionless ; here it is normalized to the 
sampling frequency, and so is given by W- 2flf I f s where f s is the 
sampling frequency. Thus, a sinusoid sampled at time t, the 
sample being the n th sample in a series of samples, can be 
expressed as either exp(jwt) or exp(jWn) . 

The data bit duration is usually several code periods 
(usually 2 0 code periods in the case of a GPS signal) . Without 
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loss of generality, the common scaling factor Ae J<p is also 
ignored since it does not affect the system design. Thus we 
assume here a normalized received signal, given by, 

x(n) = r{(n-n c )modNJe jwi "-"< ) +x noise {n) - 

Terminology 

Coherent processing: a process in which consecutive signal 
data points are summed after possibly being multiplied by 
factors that may depend on frequency; implemented as e.g. by an 
integrate and dump module followed by a filter, or by a lowpass 
filter followed by a downsampler. 

Filter (for filtering a signal) : typically a finite impulse 
response (FIR) filter, i.e. a convolution of a signal with some 
predetermined short signal (the filter coefficients) , which can 
be implemented as, 

M 

y(n) = ^h(k)x(n-k), 

where is the signal being filtered, /*(..) is a filter 

coefficient, and y(..) is the filtered signal; but could also be 
an infinite impulse response (IIR) filter (implemented e.g. 
using constant coefficient difference equations) . 

Noncoherent processing: a process in which the magnitudes 
or squares of the magnitudes of consecutive data values are 
summed . 

The Invention in General 

Referring now to Fig. 1A, in a spread spectrum receiver 
according to the invention, an initial acquisition stage 15 
analyzes a received spread spectrum signal x(t) and then 
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provides a coarse estimate w 0 of the carrier frequency and a 
value for the phase to be used for the replica code to wipe off 
(remove) the code component of the received spread spectrum 
signal. The value for the coarse estimate of the carrier 
frequency and the value for the code phase are provided to a 
fine acquisition module 16, which is also provided with the 
received spread spectrum signal x(t) . The fine acquisition 
module then determines a correction Sw 0 to the coarse estimate 
w 0 of the carrier frequency. 

Referring now to Fig. IB, the fine acquisition module 16 is 
shown as including a correlator 12 and a fine acquisition kernel 
11. The correlator 12 is typically a conventional correlator, 
of the type that would be used in a tracking module. Therefore, 
in a GPS application, the correlator would perform correlations 
of the received spread spectrum signal x(t) over one code period 

(1023 chips) . Using the coarse estimate w 0 of the carrier 

frequency and a value for the phase to be used for the replica 
code, the correlator wipes off the code component of the 
received signal x(t) and also performs a partial demodulation of 

the received signal x{t) , based on the coarse estimate w 0 of the 
carrier frequency, so as to partially remove the carrier. The 
output of the correlator 12 is therefore the navigation data 
corrupted by noise and modulated by the remaining frequency 
shift. In a typical GPS application, the bandwidth of the 
resulting signal is 1 kHz, (i.e. the correlator output will be 
attenuated for a remnant frequency modulation more than 500Hz 
from the estimate of the carrier frequency determined in the 
initial acquisition stage) . 
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The method presented here analyzes the output of the 
correlator 12 at each frequency in a predetermined set of 
frequencies to be tried as guesses at the remaining frequency 
shift. The method further demodulates the signal x(t) using 
each guess (trial value) at the carrier frequency, enhances the 
signal processed according to the guess, and uses the enhanced 
signal in a detection unit that decides whether the guess is 
approximately correct . 

1st Embodiment 

Referring now in particular to Fig. 2A, in a first 
embodiment, the fine acquisition kernel 11 is shown receiving as 
an input the approximately demodulated and code wiped received 
signal y 0 (t) , which is immediately mixed by a mixer 11a with a 
set of sinusoids, each based on a different trial value 8w 0i for 
a fine correction Sw 0 , providing the mixed samples, 

A (0 = y 0 (0exp(-y(5w 0/ 0 . (1) 

The mixed signals y u (f) (there being one for each different 
trial correction Sw 0i ) are provided to a first coherent 

processing module lib, which is implemented as either a low pass 
filter followed by a downsampler, or by an integrate and dump 
module followed by a filter. The first coherent processing 
module lib includes a (coherent) summation over a number L of 
mixed samples, and so produces one output for every L samples. 
The summation can be expressed as, 

W»)=X*/( z -» + *) (2) 

where L is the coherent integration length, i.e. the number of 
samples (mixed by mixer 11a) included in the integration 
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(summation), and n indicates the n th output, and so plays the 
role of time, but according to a different time scale than the 
corresponding index Ln+k used to indicate the different inputs 
(mixed samples) y u {L-n + k) ; the index n is incremented only once 
for every L increments of the index L-n + k . 

The first coherent processing module lib then filters the 
summation result y 2iint (n) . The filtering is performed to increase 

the coherence length (the number of samples over which the 
effect of the data modulation is negligible) without 
downsampling and so keep the signal oversampled (relative to the 
sampling that would be required for the received signal) in 
anticipation of the code wipeoff, which, as explained below, has 
the effect of doubling the remnant carrier frequency. (In other 
words, the filter further reduces the signal bandwidth, keeping 
the signal oversampled.) Although the invention does not depend 
an any particular implementation of a filter, and indeed the 
filtering here can be either a finite impulse response (FIR) 
filtering or infinite impulse response (IIR) filtering, for 
purposes of illustration, the invention is here described as 
using FIR filtering. With such a filter, the first coherent 
processing module lib provides the final output, 

y* (") = Z h 2 ( m )y2iM (n-m) = Y J h 2 (/w)£ y u {L-{n-m) + l) , ( 3 ) 

m=0 m=0 7=0 

where h 2 (..) are the filter coefficients. 

The output of the first coherent processing module lib is 
provided to a data wipeoff module 11c, which, in the preferred 
embodiment, performs essentially a complex squaring (as opposed 
to forming the product of a quantity with its complex 
conjugate) , so that in such an implementation, the output of the 
data wipeoff module can be represented as, 
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The output j/ 3| . of the data wipeoff module 11c is provided 

lid 

to a second coherent processing modul^<2"2^r, which, in the 
preferred embodiment, operates on y 3i in the same way as the 

first coherent processing module operates on its input y u , thus 
providing as its output, 



y«(n) = Z K (m)Zy^ <n-m) + l) ( 5 ) 

m=0 /=0 



where A 4 (..) are the filter coefficients, and where L 4 is the 
coherent integration length. In some embodiments, however, the 
second coherent processing is simply a coherent accumulation. 

The output y Ai of the second coherent processing module lid 
is, optionally, provided to a noncoherent processing module lie, 
which has the effect of further elevating over any noise the 
sought after signal indicating the best correction dw oi to the 

coarse estimate of the carrier frequency provided by the 
correlator 12. A noncoherent processing is, in the preferred 
embodiment, implemented as any simply noncoherent 
accumulation/summation, i.e. an accumulation (summation) that 
ignores the sign or phase of the terms being summed. Thus, 
possible expressions for the output of the noncoherent 
processing module lie include, 

ysM)-Y\y,i{L nc -n+l)\ (6) 



and 
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y^n) = f\y^L n c-n+i)\ 2 , 

1=0 



(6') 



where L nc is the coherent integration length, and where \x\ 

represents the operation of taking the magnitude of the quantity 
x , i.e. e.g. simply its absolute value if x is purely real, but 
the square root of the product of x and its complex conjugate 
if x is complex. 

The output of the noncoherent processing module lie, if 
included in an implementation of the invention, or of the second 
coherent processing module lid if noncoherent processing is not 
performed in the implementation, is provided to a detector llf . 
Whether because an implementation of the invention includes 
parallel processing for each channel / (corresponding to a 
particular trial value for a correction to the carrier 
frequency) or because the implementation processes the signal 
for each channel i in turn, the detector module llf ultimately 
receives an input corresponding to each channel. It then 
compares all of these channel inputs, each of which is the 
result of processing the signal y 0 (t) mixed with the trial value 

for the channel in such a way as to produce an input to the 
detector that is larger the closer the trial value is to the 
true correction to the carrier frequency estimated by the 
correlator 12 . In a simple implementation of the detector 
module llf, in an implementation including the noncoherent 
processing module lie, if the total number of channels is N ch so 
that the total number of trial values on each side of a null 
correction is N f =(N ch -l)/2, and A/ is the frequency interval 

separating adjacent trial values, the output of the detector 
module llf is the best estimate of the correction to the carrier 



frequency expressed as a (non-normalized) angular frequency 
increment Sw 0 given by, 

Sw 0 = 2ndf 0 = 2^{argmax(^ 50J >; 51 ,...,^ 5 ^. 1 ) - N f }-Af , (7) 
in which argmax(y 50 , j^,,...,^^,) provides the channel number for 

which the result of the processed signal is the largest, and in 
which Sf Q is the corresponding best correction to the carrier 

frequency f 0 (i.e. cycles per second). 

Referring now in particular to Fig. 2B, an implementation 
of the first . embodiment of the present invention is shown in 
which parallel processing paths are used for different estimates 
of the fine frequency correction, the embodiment otherwise being 
the same as that illustrated in Fig. 2A. The number of different 
parallel processing paths, or grid steps, is determined by the 
required frequency resolution of the acquisition system and the 
resolution provided by the correlator 12 (Fig. 1) . In the 
implementation illustrated in Fig. 2B, the first coherent 
processing module lib of Fig. 2A is implemented as a lowpass 
filter module followed by a downsample module. Thus, after the 
demodulation (by the mixer 11a in a series implementation or the 
corresponding set of mixers 11a in a parallel implementation) in 
each channel, i.e. for each frequency shift guess correction 
w 0i , the bandwidth is reduced by means of the lowpass filter 

module and then downsampled by the downsample module to remove 
redundant data from each channel . 

Referring now to Fig. 2C, in another implementation of the 
first coherent processing module lib, an integrate and dump 
module can be used, followed by a filter. 

A variation on the theme illustrated by Figs. 2B-2C is 
presented in Fig. 2D, where a DFT module is used to demodulate 
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the signal with zero-padding on the input to the DFT module and 
with output pruning (i.e. taking a subset of the output values 
for further processing) . A DFT matrix is the collection of 
sinusoidal signals of a certain length (DFT size) , each signal 
appearing as the row of the DFT matrix. In performing a DFT 
analysis of a signal, the signal is represented as a vector, the 
components of which are sample values of the signal. The result 
of the DFT analysis is the matrix multiplication of the signal 
vector. Such a multiplication involves both an elementwise 
multiplication producing a set of products, and then a summing 
of the products. These operations are similar to the 
multiplying by a sinusoid and the coherent processing of the 
invention, and so in some embodiments of the invention, DFT is 
used to implement the multiplying by sinusoids and the coherent 
processing operations. 

Complex squaring or other similar operations for performing 
navigation data wipeoff 

Because of navigation data modulation of the received 
signal, for a GPS signal the bandwidth of the correlator 12 
(Fig. 1) cannot be reduced to achieve resolution finer than 50 
Hz (corresponding to a period of 20 msec) , and the integration 
length (for coherent integration) is therefore limited to 20 
msec duration. The correlator output is, for typical GPS 
applications, at a rate of 1 msec, and the coherent integration 
length is therefore limited to 20 samples (correlator outputs) . 
Such a coherent integration length provides 50Hz resolution. As 
explained above, to obtain finer resolution, the signal samples 
are complex squared so as to wipe off the navigation data, as 
indicated by the data wipe off module 11c (Figs. 2A-2C) . But 
there are also similar operations that can be used, such as 
neighbor multiplication. 
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The terminology complex squaring is intended to convey that 
each sample is multiplied by itself (not by its complex 
conjugate) i.e. given a sampled signal x\,x 2 ,...,x h x M ,...,x N , complex 

squaring provides the signal x } x } ,x 1 x 1 ^...,x i x t ,^ M x M ,...,x N x N . The 

terminology neighbor multiplied is intended to indicate 
multiplying each sample in a sampled signal by a neighbor. More 
specifically, given a sampled signal x\,x 2 jc /+ i , neighbor 

multiplication provides ^ x 2 ,x 2 x 3 j...,^- ^p^+i^v^^-i^iv ' • Note that 
a GPS signal in a receiver is typically implemented as a complex 
number by providing by an in-phase and quadrature -phase 
component, and treating one component as real and the other as 
imaginary. 

Basis for the Invention 

The squaring doubles the modulation frequency causing 
aliasing effects. Thus the signal before the squaring device 
should be oversampled (by a factor of two) to prevent the 
frequency estimation ambiguity that could result because of 
aliasing. As already indicated, in a simple implementation such 
as that illustrated in Fig. 2C, the oversampling can be provided 
by an integrate and dump procedure followed by the trivial 
filter (i.e. the filter with coefficients [1 1] or in other 
words the filter of size 2 which has coefficients 1.0) or other 
type of filter. 

If the output of the correlator 12 is 



with d(n)-±\ denoting the data bit modulation, then after a 
coherent integration of length M, the signal will be 



y 0 (n) = A(w)d(n)e 
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y 0 '(t) = aMA{w)d(t)e J¥Vl = A\w)d(t)e 



where a is the degradation factor caused by the data 
modulation. Complex squaring wipes off the data bit, yielding, 



and so also doubles the frequency of the remaining sinusoidal 
modulation. If coherent integration is applied after squaring, 
then due to the modulation, the value resulting from the complex 
squaring will vanish if the frequency guess is wrong (w^O) but 
will not vanish if the frequency guess is correct (w = 0) , and 
may be considered to be amplified. As mentioned above, a 
noncoherent accumulation can also be used to enhance the effect, 
i.e. to further distinguish the value corresponding to the 
correct guess from the values corresponding to incorrect 
guesses. Thus, again referring to Fig. 2A-C, the invention 
includes the following steps for each channel: 

A mixing, performed by a mixer 11a, of the output of the 
correlator 11 with a sinusoid at the fine frequency shift 
correction corresponding to the channel; 

A first coherent processing performed by a first coherent 
processing module lib, i.e. lowpass filtering followed by 
downsampling (which can be implemented as an integrate and 
dump module followed by a filter) ; 

A data wipeoff performed by a data wipeoff module 11c, i.e. 
complex squaring, neighbor multiplication or other procedure 
wiping off the data modulation; 

A second coherent processing performed by a second coherent 
processing module lid; and 



j/ 0 "(0 = A a (w)d(t)e Jwt d(t)e jwt = A a (w)e 2jwt 
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Noncoherent processing (optional) performed by a 
noncoherent processing module lie, i.e. a combining of the 
magnitude of the results of the above steps over several 
periods of the carrier frequency. 

These channel steps (i.e. steps for each frequency shift guess 
correction w 0i ) are followed by a detection step performed by a 
detector module llf in which the most likely correct estimate of 
the frequency shift is made by examining the results obtained 
for each channel via the channel steps. 

The first coherent processing module lib (Fig. 2A) helps to 
_ increase the carrier to noise (C/N) ratio in a way similar to 

*n what is done by the correlator 12. Even in the presence of data 

i?! modulation, the first coherent processing module lib enhances 

N the signal. The low-pass filtering of the first coherent 

processing module lib increases the coherence length (i.e. 
ff 1 number of samples per data bit) without downsampling and keeps 

the signal oversampled. 

LI The next step, performed by the data wipeoff module 11c, 

+; (i.e. the step of e.g. complex squaring) doubles the frequency 

q of the signal by a factor of two with some loss in gain. 

Performing the step of data wipeoff by complex squaring of the 
signal input to the step preserves the phase information about 
the carrier frequency shift conveyed by the signal input to the 
step . 

The data wipeoff step is followed by a second coherent 
step, performed by the second coherent processing module lid, to 
narrow the bandwidth (frequency of the signal) down to a 
resolution of 1 Hz resolution (in case of a GPS signal) . 

What is being done here, in performing a first coherent 
processing (including a filtering) , squaring, and then a second 



-22- 



• 



coherent processing, is to first increase the signal to noise 
ratio of the correlator output mixed with a sinusoid at a trial 
frequency correction (since the signal strength increases in 
direct proportion to the integration length of the coherent 
integration, whereas the noise increases only as the square root 
of the integration length) , next to wipeof f from the signal the 
navigation data by squaring or other similar operation (which 
increases both the signal and the noise by the same factor) , and 
then to further increase the resulting signal over the noise 
using a second coherent processing. Each of the coherent 
processing operations can be viewed as performing a correlation 
of a sinusoid (at the trial frequency) with the output of the 
correlator 12 from the initial acquisition stage. Thus, the 
further correction to the carrier frequency provided by the 
invention is the frequency of the sinusoid having the largest 
correlation with the output of the correlator, and it is only by 
removing the navigation data modulation of the output of the 
correlator 12 that it is possible to distinguish between 
correlations of the output of the correlator 12 with sinusoids 
at frequencies separated by less than the frequency of the 
navigation data modulation. 

2nd Embodiment 

Referring now to Figs. 3A-C, in other embodiments of a fine 
acquisition kernel according to the invention, instead of 
demodulating the output of the correlator, as is done in the 
embodiments illustrated in Figs. 2A-D, the first coherent 
processing is performed directly on the output of the 
correlator, and the result is then complex squared (or a similar 
data wipeoff operation is performed) by the data wipeoff module 
and then analyzed for spectral content, invoking any of the 
methods of known in the art. Decoupling the coherent processing 
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from the analyzed frequency range (up to 1kHz) allows having 
several channels, each of the form illustrated in either Fig 3A, 
3B or 3C. The purpose of the first coherent stage is to reduce 
further the bandwidth by a factor of 2-L compared to the 
correlator output, and also to oversample the signal. This 
bandwidth reduction and over sampling can be accomplished by for 
example integrating over L samples and then filtering using a 
half-band filter, i.e. e.g. a filter with coefficients 
h = [1.0,1.0] . The basic idea in the approaches illustrated in 
Figs. 3A-C is to apply the first coherent processing step for 
different fine frequency estimates, square the result (to 
perform data wipeof f ) , and then analyze the frequency content 
using one or another spectral analysis tool (such as correlation 
of the result with sinusoids) . Using a single first coherent 
processing module for a group of frequencies reduces the 
complexity of the fine acquisition twofold because demodulations 
are applied at a reduced rate, i.e. the complexity is reduced 
not simply because of not providing parallel paths (one for each 
channel) until after the data wipeof f module, but also because 
demodulations are applied at a reduced rate, i.e. after the 
first coherent processing, which produces outputs at a rate that 
is down by a factor of the integration length compared to the 
input to the coherent processing. (For example, if the coherent 
processing involves an integration over 10 samples, then there 
is one output for every ten inputs.) 

For example, referring to Fig. 3A, assume that it is 
decided to apply the fine acquisition process to the range [- 
100, 100] Hz, and that the output of the correlator 12 includes 
frequencies in the range [-500, 500] Hz (i.e. including both the 
navigation data modulation and the residual carrier frequency) . 
If the first coherent processing module lib integrates 20 



consecutive samples, frequencies in the range [-25, 25] Hz will be 
filtered out by the first coherent processing process (since 
there are 20 such 50Hz intervals in the 1000Hz interval) . Each 
such filtering produces a 50Hz bandwidth channel, and four such 
channels are included in the 200Hz-wide range of [ - 100 , 100] Hz . 
Therefore, after the correlator 12, four fine acquisition 
kernels 11 as in Fig. 3A are used. The spectral estimation 
module lid' ' is used to pick out the remnant carrier modulation 
amidst noise, providing a correction to the coarse estimate 
(provided by the correlator 12) with an uncertainty down to as 
little as 1 Hz. 

The implementation illustrated in Fig. 3B, although fast, 
has as a drawback that when the integrate and dump unit is used 
for the first coherent processing, the frequency response is 
degraded for the frequencies close to the edges of the analyzed 
frequency interval because the first coherent processing module 
lib is not an ideal filter. The degradation is enhanced by the 
data wipeoff module. 

3rd Embodiment 

Referring now to Fig. 4, in another, 3rd embodiment of the 
invention, one similar to the second embodiment as indicated in 
Figs. 3A-C, a common first coherent processing module 41 is used 
to reduce the bandwidth of the signal after correlation. For 
example, if the frequency resolution of the initial acquisition 
is 100Hz, then the bandwidth of the signal after correlation is 
reduced by filtering the signal to 100Hz. Then the signal is 
demodulated as in the first embodiment (see Fig. 2A) . Thus, in 
this third embodiment, the fine acquisition kernel 11 (Fig. 1) 
will have the following modules: 



a common, coherent processing module 41a, i.e. processing the 
output of the correlator 12 (Fig. 1) before creating separate 
channels (for each trial fine frequency shift correction) ; 
and 

a filter 41b (in this case not included in the coherent 
processing module, as it is in the case of the coherent 
processing modules for the first and second embodiments) ; 

then for each channel : 

a first channel coherent processing module 42a (e.g. an 
integrate & dump module) ; 

a filter 42b; 

a data wipeoff module 42c (e.g. for performing complex 
squaring, neighbor multiplication or other similar 
procedure) ; 

a second channel coherent processing module 42d (e.g. an 
integrate & dump module) ; and 

optionally, a noncoherent processing module 42e; 

and finally, a detector 43 that scans the output of all 
individual channel computations. 

In such an embodiment, the overall coherent length, i.e. the sum 
of the coherent integration length of the common coherent 
integration 41a (before splitting into separate channels) and 
the coherent integration length of the first coherent channel 
integration 42a, should not be too large, otherwise the 
degradation due to data modulation will be high. 



The value of oversampling 

The concept of coherent processing as lowpass filtering 
plus downsampling easily generalizes to filters of arbitrary 



length. After complex squaring (or similar operation) by the 
data wipeoff module, the bandwidth of the signal is expanded by 
a factor of two. Aliasing effects then occur for a critically 
sampled signal, but can be prevented by oversampling the signal 
before the data wipeoff module using complex squaring. 

Oversampling compensates for the doubling of the bandwidth 
(resulting in aliasing) caused by complex squaring (or other, 
similar operation) of a critically sampled signal. Oversampling 
(by an appropriate factor) can provide a normalized bandwidth 
(i.e. a bandwidth divided by sampling frequency) that is 
approximately one half that for a critically sampled signal . 
Thus, the Fourier transform of the complex square (or similar 
operation) of a critically sampled signal, because of the 
doubling of the bandwidth of the signal by the complex squaring, 
will show aliasing, but the oversampled spectrum of the complex 
squared signal will not. 

Oversampling is of use not only in avoiding the aliasing 
effects that appear due to the bandwidth expansion after complex 
squaring (or other, similar data wipeoff operation) , but also 
noise amplification effects due to noise aliasing caused by the 
complex squaring. 

Comparison with zero-padding techniques 

As described above, in zero-padding for spectral analysis, 
the data of a short signal segment are padded by adding a 
segment of zeroes before applying a DFT or a fast Fourier 
transform (FFT) . The inventors have observed that although zero 
padding allows analyzing a signal at smaller frequency steps, 
and so can be considered as a tool to achieve the same ends as 
the invention, the maximum possible frequency resolution is, 
however, less than what can be achieved using the complex 




squaring of the invention even when combining zero-padding with 
noncoherent accumulation. The response of a system using zero- 
padding includes a main lobe having a width that is typically 
tens of times the width of the main lobe in the response of a 
system using the data wipeoff techniques of the invention. As 
mentioned above, the present invention reduces to zero-padding 
in an implementation that does not include a second coherent 
processing and complex squaring. 

Using different starting points 

In the above description, the method for fine acquisition 
according to the invention is used directly at the output of the 
correlator 12 (Fig. IB) . In a GPS application, the correlator 
output provides samples of the signal at 1ms intervals, and 
there are 20 such samples per data bit. It is assumed that the 
timing of the bit transitions (i.e. when they occur) is unknown, 
and the inventors have verified the method of the invention 
based on such an assumption. Thus, according to the preferred 
embodiment of the invention, the processing after the correlator 
starts somewhere in the middle of a data (navigation) bit (i.e. 
somewhere in the middle of a 2 0 -sample interval) at one or 
another sample. 

However, in another embodiment of the invention, the method 
starts from two or several different samples. In such an 
embodiment, the method first starts as in the preferred 
embodiment, at one or another sample, and proceeds as described 
above. But then the some predetermined number of samples (e.g. 
10 samples) are skipped, and the method is restarted at a new 
starting point and again proceeds as described above. The 
procedure is performed for each of selected samples. In such an 
embodiment, therefore, the "phase" of the starting point with 
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respect to the navigation data bit edge is changed. Such an 
embodiment might occasionally provide increased performance. In 
a straightforward implementation of such an embodiment, the 
processing for the two or several starting points is performed 
seriatim. In an alternative implementation, one designed to 
reduce acquisition time, the processing for the two starting 
points is performed in parallel. 

SCOPE OF THE INVENTION 

It is to be understood that the above-described 
arrangements are only illustrative of the application of the 
principles of the present invention. In particular, some of the 
computation required by an acquisition system according to the 
present invention can be performed outside of the device 
receiving the spread spectrum signal, such as by a component of 
a network using the spread spectrum signaling. Moreover, fine 
acquisition according to the invention can be used independently 
of an initial, coarse acquisition; it is useful in any situation 
where a signal includes a carrier component that must be 
estimated, but does not include a code component (i.e. in case 
of a code division multiple access signal, the code is wiped off 
the signal before the signal is processed according to the 
invention) . In addition, other parallel and sequential 
architectures are comprehended besides what are expressly 
indicated herein, and numerous other modifications and 
alternative arrangements may be devised by those skilled in the 
art without departing from the spirit and scope of the present 
invention, and the appended claims are intended to cover such 
modifications and arrangements. 
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